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ABSTRACT 


Exact  solutions  are  presented  to  a  number  of  small  perturbation 
magnetohydrodynamic  flow  problems.  The  conditions  under  which  the 
solutions  are  obtained  are  as  follows: 

1.  The  flow  is  two-dimensional,  and  is  only  slightly  perturbed 
from  a  uniform  flow. 

2.  The  magnetic  field  vector  is  also  two-dimensional  and  lies 
in  the  plane  of  the  flow. 

3.  The  distortion  of  the  applied  field  by  the  induced  currents 
is  neglected. 

4.  Physical  boundaries  on  the  flow  are  one  or  two  infinite 

plates  parallel  to  the  flow  direction. 

5.  The  conductivity  of  the  fluid  is  a  scalar  quantity,  but  may 
vary  with  position. 

With  these  assumptions,  the  perturbations  to  the  flow  are  calculated 
for  various  magnetic  fields  (chiefly  those  due  to  a  current  flowing  in  a 
single  wire,  and  a  linear  dipole)  for  incompressible,  subsonic  and  super¬ 
sonic  free  stream  speeds.  Calculations  of  the  pressure  on  the  walls  and 
other  quantities  are  presented  for  illustrative  examples,  including  cases 
in  which  the  conductivity  is  not  uniform  throughout  the  flow. 
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LIST  OF  SYMBOLS 


Dimensional  quantities  may  be  taken  to  be  expressed  in  the  rationalized 
MKS  system  of  units. 

Dimensional  Non-Dimensional 


4* 

J 

current  vector 

j*z 

z- component  of  the 
current  vector 

<7* 

a1 ,  a 

conductivity 

velocity  vector 

U*,  V* 

u',  v';  u,  v 

velocity  components 

u 

free  stream  velocity 

B* 

magnetic  field  vector 

B*  .  B* 
x  y 

b  ,  b 
x  y 

components  of  the 
magnetic  field  vector 

B 

o 

reference  magnetic 
field  strength 

x*,  y* 

x,  y 

rectangular  Cartesian 
coordinates 

reference  length 

p* 

p',  p 

mass  density 

Po 

reference  mass  density 

p* 

P'>  P 

pressure 

Po 

reference  pressure 

s* 

s',  s 

entropy 

T* 

temperature 
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Dimensional 


N on  -  Dimens  ional 


Y 

S 

M 

£  *  n,  t,  £ 

z 

w 


*o 

I 


contours. 


P 

H,  X 
X  .  fi 
* 


F(x) 


gas  constant 

ratio  of  specific  heat 

interaction  parameter 

Mach  number 

dummy  variables 

complex  variable 

complex  function  of 
complex  variable 


lengths 

complex  variables 

special  function 

stream  function 

permeability  of  free 
space 

current 

complex  function  of 
real  variable 


Subscripts  Pand  c  refer  to  particular  and  complementary 


solutions. 
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1 .  Introduction 


A  number  of  interesting  steady  flow  problems  in  inviscid  magneto¬ 
hydrodynamics  may  be  treated  approximately  by  neglecting  entirely  the 
magnetic  Reynolds  number,  and  assuming  further  that  the  interaction 
parameter  (as  defined,  say,  in  Ref.  1)  is  small.  In  this  paper  solutions 
are  given  to  a  number  of  these  problems  under  the  following  additional 
restrictions:  first,  the  flow  is  two-dimensional;  second,  the  magnetic 
field  lies  in  the  plane  of  the  flow;  third,  the  flow  is  only  slightly  perturbed 
from  a  uniform  flow,  the  perturbation  being  of  the  first  order  in  the  inter¬ 
action  parameter. 

Possibly  the  simplest  flow  of  this  type  is  that  considered  by  Sherman; 
the. flow,  which  in  this  case  is  taken  to  be  incompressible,  is  through  an 
infinite  channel  of  arbitrary  height,  while  the  magnetic  field  is  due  to  a 
current  flowing  in  a  single  wire  located  outside  the  channel  and  running 
perpendicular  to  the  flow  direction.  This  problem  was  treated  numerically; 
the  resulting  pressure  gradient  on  the  surface  of  the  channel  was  used  in  a 
subsequent  paper 3  to  investigate  the  boundary  layer  over  one  wall  of  the 
channel. 


A  similar  problem  in  the  sense  of  being  a  small  perturbation  on  a 
uniform  flow  was  discussed  by  Kemp  and  Petschek,  4  but  the  geometrical 
configuration  treated  in  that  paper  was  somewhat  different. 


It  will  be  shown  in  this  paper  that  the  solution  to  the  channel  problem 
may  be  obtained  analytically;  the  case  in  which  the  channel  height  is  infinite 
has  a  simple  closed-form  solution.  In  addition,  a  method  will  be  given  for 
solving  similar  channel  flow  problems  for  any  two-dimensional  magnetic 
field  arising  from  an  arrangement  of  pole  pieces  or  current  carrying  wires 
outside  the  channel.  Finally,  solutions  are  given  to  some  linearized  sub¬ 
sonic  and  supersonic  compressible  flow  problems  with  the  same  geometry, 
and  some  incompressible  problems  in  which  the  conductivity  is  not  uniform. 


2.  General  Analysis 

In  the  absence  of  an  applied  electric  field  the  current  is  given  by 
the  expression 

B* 

(2.  1) 

Taking  Cartesian  coordinates  (x*,  y*)  in  the  plane  of  the  flow,  with 
the  x-axis  parallel  to  the  undisturbed  flow,  it  is  found  that  the  current  is 
entirely  in  the  z -direction  and  is  given  by 


=  o-  (u.*  6^  -  v  B*  ) 

The  body  force  on  the  fluid  is  given  by 

and  the  equations  of  motion  become 


(2.2) 


(2.3) 


(2.4) 


<*«**?♦ 

(2.5) 

In  addition,  the  equation  of  continuity  is 

•£*0**u 

(2.6) 

and  the  energy  equation  is 


(2.7) 

Finally,  we  assume  that  the  fluid  is  a  perfect  gas  with  constant  ratio  of 
specific  heats;  so  that 

p*.p*RT* 

(2.8) 

(2.9) 

The  following  non-dimensional  quantities  will  now  be  defined: 


w*  y*  .  j# 

,  Y*  ■£■  ,  j‘o»d?;  , 
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(2.  15) 


d»'« 


Where  S  (the  interaction  parameter)  is  defined  by 


(2.16) 

The  magnetic  field  is  assumed  to  be  unaffected  by  the  currents 
flowing  in  the  gas  (low  magnetic  Reynolds  number);  this  assumption  will 
always  break  down  far  from  the  sources  of  the  field,  but  will  be  satisfactory 
for  the  purposes  of  this  paper.  In  these  circumstances,  Maxwell's  equations 
will  require  the  magnetic  field  to  satisfy 


(2. 17) 


(2.  18) 


Subject  to  these  conditions,  however,  the  magnetic  field  will  be  a  given 
quantity  in  any  particular  problem. 


The  flows  to  be  considered  in  this  paper  will  be  only  slightly  per¬ 
turbed  from  a  uniform  parallel  flow;  perturbation  quantities  are  defined 
as  follows: 


w.'  *  l  +■  Su. 


(2.  19) 


v'*$v 


(2.20) 
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(2.21) 
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The  last  equation  may  be  integrated  at  once  to  give 


e 


(2.29) 

Equations  (2.  24)  -  (2.  27)  and  (2.  29)  present  the  problem  under  considera¬ 
tion  in  a  form  suitable  for  analysis,  but  before  considering  any  particular 
problem,  some  general  conclusions  may  be  drawn.  First,  Eq.  (2.  27)  may 
be  integrated  at  once  to  give 


su,y)  =  rMa-/XbviCy)d$ 

-80  7 


(2.  30) 


Second,  the  linearized  momentum  equations,  (2.  24)  and  (2.  25),  now 
describe  a  flow  moving  under  a  specified  non-conservative  body  force;  this 
force  is  calculated  directly  from  the  given  field  quantities.  Third,  the 
x-momentum  Eq.  (2.  24)  may  be  integrated  to  give 


by*<g,y)d£ 

*'-80 

(2.31) 

Next,  the  terms  in  u  and  p  may  be  eliminated  from  Eq,  (2.  26)  by 
using  Eqs.  (2.  24),  (2.  27)  and  (2.  29)  to  give: 


+(mz-  0  t + or- 0  M*] byx 

(2.32) 

Now,  if  Eqs.  (2.  25)  and  (2.  32)  can  be  solved  for  v  and  p,  u  may  be 
found  from  Eq.  (2.  31),  and  p  from  Eqs.  (2,  29)  and  (2.  30),  thus  completing 
the  solution  of  the  problem. 

3.  Incompressible  Problems 


In  the  incompressible  case  Eqs.  (2.  25)  and  (2,  32)  may  be  written 

37*  §y  = 

(3.  1) 
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4*-  ~  *  Li 

ax 


which  equation  together  with  Eq.  (2.  31): 


(3.2) 


/* 

U.s-p-1  by^yjdx 

•OO 

define  the  problem. 


(3.3) 


It  is  a  remarkable  fact  that  a  particular  solution  of  Eqs.  (3.  1)  and 
(3.  2)  may  be  found  for  an  arbitrary  magnetic  field.  This  solution  is: 


vp.iby/\«,vW« 


(3.4) 


“<© 


(3.5) 


The  integrals  in  these  expressions  will  be  convergent  in  all  normal 
cases.  The  complete  solution  to  Eqs.  (3.  1)  and  (3.  2)  can  now  be  written 


vp  +  vc 


F>=Pp+Fc 


where  vc  and  pc  satisfy 


(3.6) 


a*  dr  0 


^-4&=0 

ay  dx 


(3.7) 
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Clearly,  Eq.  (3.  7)  may  be  satisfied  by  writing 


p  +i  vc  =w6c+iy)  =  wfe) 

C  (3.8) 


where  w  (z)  is  an  analytic  function  of  the  complex  variable  z. 

The  boundary  condition  appropriate  to  any  particular  problem  will 
require  the  normal  component  of  the  velocity  to  vanish  at  a  solid  surface; 
in  this  paper  we  will  consider  only  infinite  channels  defined  by  y  =  0  and 
y  =  H,  say,  together  with  the  limiting  case  H-*  °° .  On  these  surfaces  the 
boundary  condition  then  reduces  to 


(3.9) 


Single  Wire  Case 

Turning  now  to  particular  cases,  the  simplest  magnetic  field  is 
that  due  to  a  current  I  in  an  infinite  wire  through  the  point  (0,  -1)  with  the 
lower  wall  of  the  channel  at  y  =  0.  The  physical  distance  of  the_wi^e  from 
the  wall  is  y  ,  and  B  may  be  taken  to  be  the  field  at  the  origin,  *A°  /Zny  . 
Then  °  ° 


2. 


xi+<Y+0*' 


and 


V 


P  = 


4*W  vrr 

X*  -t  (Y+  |V- 


=  j[x-(y+i)t*r,->  ^7-1 
rP  x*  ♦  Cv+  0*  " 


(3.10) 


(3. 11) 


(3.  12) 

In  deriving  these  expressions  from  Eqs.  (3.  4)  and  (3.  5)  all  the 
integrals  were  taken  from  the  lower  limits  C,  =  0  instead  of  £,  =  -°o.  It  may 
be  verified  that  in  this  particular  case  Eqs.  (3.  11)  and  (3.  12)  are  still  par¬ 
ticular  solutions  to  Eqs.  (3.  1)  and  (3.2).  A  slight  simplification  of  the 
algebra  results  from  this  unimportant  modification. 

Considering  first  the  case  of  infinite  channel  height,  w  (z)  must  be 
analytic  for  >  0,  and  satisfy: 
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on  y  =  0 


*0  w(z)  = 


X  tdw"'x 
X*  +1 


(3.13) 

The  solution  may  be  obtained  by  standard  methods  and  is  easily  shown  to  be 


w<2)= 


>  -  iz 


(3.  14) 

The  complete  flow  field  is  now  easily  established,  but  it  will  be 
sufficient  here  to  derive  the  more  important  features  of  this  flow.  The 
pressure  on  the  wall  is  given  by 


PCX,0)» 


X-Wx  -fex  1„  ^ 


(3.15) 

In  addition,  we  find 


i  L»Y 

v(o.v)=  — 

2.  (y*-i) 

(3.  16) 

and 


Zu-(o,y)=  w.(oo  y)=-^_ 


(3.17) 


This  profile  is,  of  course,  rotational,  and  represents  the  amount  of 
vorticity  produced  by  the  non-conservative  force  field.  Finally,  the  pres¬ 
sure  gradient  on  the  wall  is  found  from  Eq.  (3.  15)  to  be 


.  (**+■)* 


(3.18) 


At  the  origin  this  gives-^-(0,  0)  =  ln'TY.  Equations  (3.  15)  and  (3.  18)  are 
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shown  in  Figs.  1  and  2,  together  with  material  to  be  described  later. 

Single  Wire  and  Channel 

Next,  consider  an  infinitely  long  channel  with  walls  at  y  =  0  and 
y  =  H,  the  magnetic  field  being  the  same  as  in  the  previous  case.  vp  and  pp 
are  still  given  by  Eqs.  (3.  11)  and  (3.  12),  but,  in  addition  to  Eq.  (3. 13), 
w  (z)  must  satisfy 


w  s 


"^FTThTo27 


on  y=H 


(3.  19) 


Defining 


I  -  \z 
Z  H 


t  +iz 
1  H 


(3.20) 


w  (z)  may  be  shown  to  be  given  by: 


w(2)s  Itkia+il^oth  -  cothi-irA-  2-]- 


_ L 


?Hi  n=  i 


In  (w-i-m) 

(n-XKn-i-X) 


Cn-MKn-t-p) 


(3.21) 

Interesting  deductions  that  can  be  derived  from  Eq.  (3.  21)  without  summing 
the  series  are: 


CoW)  —  u.  (-o.Y)  =  ^ 


lw(H-H) 

M 


2P(0,V)=  p(c©,y)  =  -TrhjH  +  O 


(3.  22) 


(3.23) 
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This  expression  for  the  pressure  drop  along  the  channel  is  plotted 
in  Fig.  3.  For  finite  channel  flows  of  this  type,  the  drag  is  just  the  product 
of  the  channel  height  and  the  pressure  drop.  However,  it  is  not  necessary 
to  calculate  the  flow  field  to  determine  the  drag;  to  this  order  it  is  only 
necessary  to  integrate  the  magnetic  force  through  the  field.  When  H-*-  °o 
the  drag  is  logarithmically  infinite,  a  result  which  is  due  to  the  slow  way 
in  which  the  field  decreases  at  large  distances  from  the  wire.  The  pressure 
on  the  lower  wall  (y  =  0)  is  given  by 


Pfr,o)s-:ZLkffl±i3 


shx^ 


^(x-tan'x) 

x’-TT 


i*  In 


x*±i 


4*""  )*+<*>*! 


(3.24) 

This  expression  is  shown  for  various  values  of  H  in  Fig.  1,  and  its 
derivative  in  Fig.  2.  It  can  be  shown  that  for  any  finite  H  the  pressure 
tends  to  its  asymptotic  value  from  above,  while  for  infinite  H  the  reverse 
is  true--an  interesting  case  of  non-uniform  convergence.  The  curves  for 
H  =  0  correspond  to  the  assumption  v  =  0  (from  continuity)  and  come  from 
integrating  Eq.  (3. .3)  on  this  basis. 

The  results  obtained  for  the  pressure  gradient  appear  to  agree  quite 
closely  with  those  of  Sherman;3  the  actual  pressures  differ  by  a  nearly 
constant  amount;  this  is  probably  due  to  the  fact  that  the  numerical  inte¬ 
gration  was  commenced  by  setting  p  =  0  at  some  finite  upstream  distance 
and  is  of  no  importance.  The  pressure  gradient  is  unfavorable  for  a 
boundary  layer  on  the  wall  over  a  distance  of  the  order  of  unity  near  the 
wire.  This  unfavorable  region  shrinks  with  the  channel  height. 

Two  Wires 


As  an  example  of  a  slightly  more  complicated  magnetic  field,  as 
well  as  a  more  realistic  case  in  the  sense  that  the  net  current  giving  rise 
to  the  magnetic  field  is  zero,  consider  two  wires  located  at  the  points  (X,  -1), 
(-X,  -1)  carrying  equal  currents  I  in  opposite  directions.  The  distance  of 
either  wire  from  the  wall  y  =  0  is  yQ  and  the  distance  between  the  wires  is 
2Xy0.  B0  will  again  be  taken  to  be  -  p.Ql/27ryo,  so  that  the  non-dimensional 


field  components  are 


k  -  v»* 

c^~Cx+-)6i-+(yti)i  “ 

(3.25) 

The  particular  solution  in  this  case  is: 


<y»o 


<x-X> 

WZlf+O* 


so  that  for  the  single  wall  case  w  must  satisfy 


(3.27) 


^  =  i  VX>  W<x*>)[ -  -g***]  on  y=o 

The  solution  to  this  problem  is 


(3.28) 


(3.29) 
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Hence 


-2W-(0,y)aU.<oo,V)» 


-•ft-X* 

(Y  +  OLX^+CY+t)*^ 


p(o,y)=o 


(3.30) 

(3.31) 


P(,X,0?= 

(x+X>*+ »  (x-X>*-  +  i  iFX-l"  (x+X)*  +  » 


(3.32) 

The  drag  in  this  case  is  finite;  it  can  be  found  from  integrating 
Eq.  (3.  30),  but  it  can  also  be  found  from  integrating  the  force  in  the 
x-direction  through  the  whole  flow  field  without  calculating  the  actual 
velocity  components.  Either  way,  the  result  is-jp-  In  (1  +  X2).  The_dimen- 
sional  drag  per  unit  length  is  then  crUB02y02.  -2-  In  (1  +  X2)  =  crUp^I2  . 

A  In  (1  +  X2).  The  pressure  on  the  wall  Eq.  (3.  32)  is  shown  in  Fig.  4  for 
a  number  of  values  X;  the  pressure  gradient  is  shown  in  Fig.  5. 

For  comparative  purposes  the  product  of  the  current  in  either  wire 
times  the  distance  between  the  wires  has  been  kept  constant.  This  has  the 
effect  of  making  the  reference  field  equal  to  the  y-component  of  the  field 
at  the  origin  and  permits  transition  to  the  case  of  the  linear  dipole.  For 
large  X  the  two  wires  behave  independently  as  may  be  seen  by  comparing 
with  the  previous  case.  However,  as  X  gets  smaller  the  character  of  the 
pressure  distribution  changes;  the  pressure  gradient  becomes  unfavorable 
over  a  longer  distance,  but  has  smaller  values.  The  curves  for  X  =  0.  1  are 
indistinguishable  from  those  for  the  dipole  case  (X  =  0)  which  will  be  presented 
next. 

Linear  Dipole 

Particularly  simple  results  may  be  found  in  the  case  of  the  linear 
dipole.  If  this  is  placed  at  the  point  (0,  -l)(distant  yQ  from  the  wall)  and 
the  magnetic  field  at  (0,  0)  is  taken  to  be  B0,  the  non-dimensional  com¬ 
ponents  of  the  magnetic  field  are: 


-13- 


u- 

**  [x^+ty+O1]*- 


L  x*-cy -H)*- 

V  [xSCy*«)*l*- 


(3.  33) 

This  case  may  also  be  derived  from  the  previous  case  by  letting  I  become 
large  and  X  become  small  in  such  a  way  that  2X1  =  27rBoyo/pQ.  The  par¬ 
ticular  solution  is: 

VP=  [x*  +  <y+0*]* 


VxtCy+if-KTX*] 
fxz+Cyt  trP 


so  that,  for  the  infinite  channel  height  case, 


n  _-Jt(XX-|) 


on  y sQ 


(3.34) 


(3.35) 


(3.36) 

Hence 

w  * -£■  (z*>3  ia-4Xz+iV? 

(3.37) 

The  complete  solution  is  then: 

-a(?+W  -—■)  | x [x1  (y*-+2y-i) t (y+ 1 )2'(yz+4Y-0] 

U"  (y+03  +  *  (y+ 


v=  ^y^-*-(y-H)(yi-,3)] 


[x1+(y+ 1  )*•]*• 


(3.  38) 


(3.39) 
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_  - j-x  [** -V  (y*+  3j  +?)+  (y+Q^y*-*  +y+s  >] 


Also,  the  stream  function  defined  by 


u. 


is  found  to  be 


.i  yCy+-o 

8  trnw 


-)-± 


XY 


(3.40) 


(3.41) 


The  pressure  gradient  on  the  wall  is 


dx 


pCx,o)  = 


_  xb-3tt4+23x Zr  £ 


% 


■0* 


(3.42) 


(3.43) 

Lastly,  the  streamline  which  is  defined  by  i|i  =  - yQrises  in  passing  the  dipole, 
the  amount  of  the  rise  being 


Ay  = 


Sir 

~8~ 


Cyo+0*- 


(3.44) 


The  pressure  and  the  pressure  gradient  are  shown  in  Figs.  4  and  5. 
Linear  Dipole  and  Channel 


If  a  finite  channel  is  introduced  having  walls  at  y  =  0  and  y  =  H,  the 

(3.45) 


conditions  on  are  Eq.  (3.  36)  and 

*0  w  =  -■^(H+0[xt-(H+i)2y£xi+  (H+t)*}2 


The  complementary  solution  is  then 
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(3.49) 


The  last  expression  gives  the  pressure  drop  along  the  channel  and  is 
shown  in  Fig.  6.  The  pressure  on  the  wall  y  =  0  is 


P(x,o)= 


Tt'CH+a.) 

ICtH+Q*- 


X 


£ 

«=» 


* 


(3.  50) 

and  is  shown  in  Fig.  7.  The  effect  of  the  channel  in  this  case,  as  in  the 
case  of  the  single  wire  is  again  to  reduce  the  region  of  unfavorable  pressure 
gradient;  eventually,  as  H-*0  the  pressure  drops  continuously  along  the 
channel. 
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4.  Subsonic  Problems 


No  solution  of  the  generality  of  Eqs.  (3.  4)  and  (3.  5)  is  known  in  the 
case  of  subsonic  compressible  flow.  It  is  not  sufficient,  for  the  inhomoge¬ 
neous  problem,  merely  to  scale  the  quantities  according  to  the  Prandtl- 
Glauert  rule.  This  is  because  the  existence  of  the  particular  solution  Eqs. 

(3.  4)  and  (3.  5)  depends  crucially  on  the  relations  Eqs.  (2.  17)  and  (2.  18) 
connecting  the  components  of  the  magnetic  field;  these  relations  are  modified 
in  the  scaling  procedure  leaving  the  problem  in  a  new  form  that  is  apparently 
no  simpler  than  the  original  one.  However,  particular  solutions  can  never¬ 
theless  be  found  for  simple  field  configurations;  once  such  a  particular 
solution  is  found,  for  some  magnetic  field,  the  homogeneous  problem 
becomes  a  straightforward  problem  in  complex  variables  whatever  the 
nature  of  the  boundaries. 

Single  Wire  Case 

First,  take  the  magnetic  field  as  due  to  a  current  in  single  wire,  the 
components  being  given  by  Eq.  (3.  10).  The  equations  for  v  and  p  are: 


dv  dp 

dx  dy 


-x£y+_1> 
[x^+Cy  +,yf 


(4.  1) 


(4.2) 


By  transforming  to  polar  coordinates  with  origin  at  the  wire,  it  can  be  seen 
that  particular  solutions  of  Eqs.  (4.  1)  and  (4.  2)  exist  in  which  v  and  p  each 
have  the  form  of  a  function  of  the  angle  divided  by  the  first  power  of  the 
distance ;ion  transforming  these  particular  solutions  back  into  Cartesians, 
one  obtains: 


.  ^-ECi+U-DM^xCtan'^rr)  . 

'P - x^TFCy-h^  *  °  [** 


+(Y+0i 


x*"Ky+i)2 


Y 

Xz+/S1(Y+0z 


(4.3) 


where 


for  subsonic  flow. 


(4.  4) 


The  complementary  solution  vc,  pc  satisfies: 


dy 


=  o 


dvc.  . 
ax  + 


-o 


(4.5) 


(4.6) 

and  these  equations  will  be  satisfied  by  an  analytic  function  w  (z)  defined 
by:  c 


wc(zsx+ip  y)  =  (Jpt+iVc 


Taking  the  case  of  a  single  wall  at  y  =  0,  the  condition  on  vc  is: 

*l  '  ’  2r[F7T 


The  solution  of  this  problem  is 


W, 


■  '  +  +  £*(z+i  Z+»eJ 


(4.7) 


(4.8) 


(4.  9) 


From  this  may  be  found  all  the  remaining  flow  quantities.  For  example: 

ifln  \'^\r0) 


(4.  10) 


Hence,  from  Eq.  (2.31) 
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Furthermore,  from  Eq.  (2.  30), 


six,o)*irMx[(|f+u«ii)- 


(4. 11) 


and  from  Eq.  (2.  29) 


eMU  —  lz  ? ln  (T^y 


+  t4h  X 


rScUi*™S£ 


(4.  12) 


(4. 13) 


These  quantities  are  shown  in  Figs.  8,  9.  10  and  11  for  various  values  of 
M,  and  for  y  =  1.4  and  y  =  1.0.  Comment  on  these  graphs  will  be  reserved 
until  after  the  discussion  of  the  supersonic  case  which  will  be  given  in  the 
next  section. 

The  channel  problem  can  also  be  solved  by  the  method  of  images, 
but  the  solution  will  not  be  presented  here. 

Linear  Dipole 

The  transformation  into  polar  coordinates  can  also  be  used  to  find 
the  particular  solution  to  the  subsonic  problem  when  the  magnetic  field  is 
due  to  a  linear  dipole,  the  components  of  the  magnetic  field  being  given  by 
Eq.  (3.  33).  The  particular  solution  is: 


PP - 4lxx_+(?i(y+0AP  '  Cx-taw  y& 
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+  3(8  -S?Z~  if +)x\v^tf+(9-  I30Z-  .24^  xVn)-  3  |  +  3  /$>*(Y+»)  ®| 

*  lUx^m.V^o*]3  **>*<»>? 


4fl- (3 pZ-  4^ V(V+D-3^(H- 3/8J*)(y+.> 


1 


xlx^iytif]  r..„.,»uti  M*/TT  .*_-'  x  V 

'  D^^Wt  l  )1  4("S-t“  vn> 


(4.  14) 


-vtit-O  M*x*(y+  0 _ 

U-txZ4  <Y+  0*)3  txi+^?-  (v+  l)*]3 


JV(  1  -/5*-)Xfc+(4- 23  0*'- S/S*)  X*(y+0* 


43(M3/S*- 4^4)xHy+iA  3^iC3+5^;iXv+-<)^ 


-*1/3*4  ;L  /«  ^(Y*  0*>  (3  '31  fi3--  <v84+  Cy+ 1 A  3^(-3-5/<J'V4/tf +)(y+  i)^J 

(4. 15) 

The  homogeneous  problem  can  again  be  solved  by  a  modified  complex 

variable  as  defined  in  Eq.  (4.  7).  Taking  the  case  of  a  single  wall  at  y=  0, 

the  condition  on  v  is: 

c 

VC  cx*0)  =~  [’  +  Cv‘  ‘>  M*]  T"  (?  -  taw"  x)- 

-  [4 ( »  -£*)**>  w- -*3**  J^+)X4+3(t- i30*-+/P)X%30\~3-S/8Z) j 


I [3^*> <» '•-itfif+lfiW 


+(3-3,l/$*'-t>/3\t>/3i)x:L+3/S2'('3~S/3l++/2i')  J 


(4.  16) 
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The  solution  to  the  homogeneous  problem  is  then 

WC  .[.♦<*-  o  M*]  I  +  5?;^ 


•(z*-3^^1h  T Tjt~ 

zz+/3^ 


t  £jr£<± ig£?1  - [j-  -  +A i >  ^  _  3  _> 

I(h^F  J  ,2-  15^  (*+'>*  <z+Oj  %  z  +  t/3 

v  z&Mh.  *±u£.  -  [3  _J _ 

*  zfi  (z+W:  Xfi  ^+'/313  J  <2-  [1-76*  2+i 

■3(Z+*Q  ■  3zz+<)iz-  t _ 3  _J _  .  3U+Z&*)  z+Zi  fi 

(*+0*  (Z+TP  /S<»-/3;*  2+1/6 +  *,<?  u7i£p 


(z+ij^  T  (z  +  i)-3 

-  (3+ a/?*)  foz^+^'zg-aff^ 
2<j  ^z+r&p 


y 

m 


(4.  17) 


and  the  pressure  on  the  wall  may  be  shown  to  be: 


M* 


pc*,o>=[i+ar-i)M2]  -jjj- 


(x**/S*)4  +  Tx*+ 


.  xx-tfl&+Sfi)  (  Y-I)x  ~3  v  2  — 8 1*JL_ 

+  Z(i +/5)*  J  t  IZ/3  ^  '  ;L^+»  (Xi+0a*  CXi+»)aJ  I2-/S 

[  Gift's  +  0<W] 


Proceeding  as  before, 


*M*,o)  =  -PCx»o)-i  (ST  +tdrf  >0- 


i  »  +1  * — 

+  X*+l  6  (**+ 1) 


(4.  18) 


(4.  19) 


(4.20) 
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(4.21) 


f>(x,o)*M*P-  1X^1  S 


These  quantities  are  shown  in  Figs.  12,  13,  14,  and  15  for  various 
values  of  M  and  for  y  =1.0,  and  y  =  1.4.  Comment  on  these  graphs  will 
also  be  reserved  until  comparison  can  be  made  with  the  supersonic  case. 

It  is  clear  that  the  algebraic  complexity  of  the  subsonic  solutions 
will  preclude  the  use  of  more  complicated  magnetic  fields.  On  the  other 
hand,  the  two  problems  which  have  been  solved  serve  to  show  the  trend  of 
the  effects  of  compressibility  at  subsonic  speeds  on  the  type  of  flow  being 
considered.  It  is  interesting  to  observe  that,  as  in  ordinary  linearized 
compressible  flow,  the  pressure  increases  with  (1  -  M^)"*'  . 

5.  Supersonic  Problems 

When  the  flow  is  supersonic,  the  problems  under  consideration  can 
be  treated  by  the  method  of  characteristics.  Introducing 


£=|£x-  x(Y+0,^“i£x*ir(Y+l) 


(5.  1) 


(5.2) 


2  2 

where,  for  supersonic  flow,  (3  =  M  -  1,  Eqs.  (2.  25)  and  (2.  32)  become: 


•^-[0P-v)  -“-pf  ^fev+0+<5M)Mz}  &  by 


(5.3) 


[0p-t-v]=  ^  ^ 


(5.4) 
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When  there  is  just  a  single  plate  at  y  =  0,  the  problem  can  be  solved 
for  arbitrary  field;  the  appropriate  characteristics  are  shown  in  Fig.  16. 

In  order  to  find  the  conditions  on  the  plate  it  is  only  necessary  to  integrate 
Eq.  (5.  3)  along  the  characteristic  marked  1;  since  v  (x,  0)  =  0, 


?(x.a)=-  I  ^ 

0© 


4r  0 It).  j  by  j^+«>  |  * 


.  r»+tv-0Mz  n-* 


V-t  -  tA 

M 


<k 


(5.5) 

In  order  to  find  conditions  anywhere  else  in  the  flow  field,  two  additional 
integrations  are  required,  one  from  infinity  along  the  characteristic  marked 
2,  and  the  other  from  the  plate  (where  v  =  0  and  p  is  given  by  Eq.  (5.  5)) 
along  the  characteristic  marked  3.  By  way  of  illustration,  the  two  cases 
considered  in  the  last  section  will  be  worked  again  for  the  supersonic  case. 


Single  Wire  Case 

For  the  single  wire,  with  the  components  of  the  magnetic  field  given 
by  Eq.  (3.  10),  the  pressure  on  the  wall  becomes: 

]  X+£ 


(5.  6) 


where  sgn  (x  +  p)  =  +1,  0  or  -1  according  as  (x  +  p)  is  positive,  zero  or 
negative.  The  rise  in  pressure  at  transonic  Mach  numbers  is  again  evident. 
Proceeding  as  before, 


U.CX, 


(5.7) 


$(x />)=  ^  y  M*  [(£+&.*)' ‘x)~  — 


(5.8) 
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-4or-OM*,(f +u«*‘x) 


(5.9) 

These  quantities  are  shown  on  Figs.  8,  9,  10  and  11  along  with  the  subsonic 
solutions.  The  following  points  should  be  noticed  in  connection  with  these 
graphs.  First,  the  curves  for  y  =  1.  0  may  be  thought  of  as  emphasizing 
the  force  aspects  of  the  flow,  while  the  curves  for  y  =  1.  4  bring  out  the 
effect  of  heating.  This  is  particularly  clear  in  Figs.  11a  and  lib  where  it 
will  be  noticed  that  the  change  in  density  is  much  greater  for  the  case  M  =  2, 
■y  =  1. 4  than  M  =  5,  y  =1.0.  This  may  be  thought  of  as  being  due  to  the 
fact  that  large  Mach  numbers  correspond  to  low  temperatures;  in  this  case 
conversion  of  a  small  fraction  of  the  flow  energy  into  heat  is  sufficient  to 
cause  a  large  rise  in  the  temperature  of  the  flow  accompanied  by  a  large 
expansion.  From  9a  and  9b  it  can  be  seen  that  the  change  in  velocity  is 
largely  independent  of  the  Mach  number,  thus  the  flow  energy  converted 
into  heat  is  roughly  constant  but  has  a  much  greater  effect  on  the  flow  at 
high  Mach  numbers.  Finally,  this  phenomena  will  set  a  limit  on  the  validity 
of  these  calculations  as  follows:  a  small  change  in  the  temperature  will 
produce  a  relatively  large  change  in  conductivity  (cr  normally  varies  as 
quite  a  high  power  of  T).  Thus,  for  high  Mach  numbers  A  T/T  cannot 
exceed  some  quite  small  value  in  order  that  the  assumption  of  constant 
conductivity  should  remain  valid. 

As  far  as  the  pressure  effects  are  concerned,  the  dominant  effect 
is  transonic;  a  large  rise  in  pressure  appears  as  M  -*•  1  for  both  the  super¬ 
sonic  and  transonic  cases;  this  is  a  characteristic  phenomenon  in  all  com¬ 
pressible  flows,  and  its  appearance  in  this  magnetohydrodynamic  case  is 
quite  interesting. 

Linear  Dipole 

For  the  linear  dipole,  the  quantities  are  given  by: 


P(X,0)  = 


(X3-  3/SXL-3K+/3) 


[i+tK-OM*] L  -‘/Sx-i.-n r  (x+/<S> c^x-o  1 
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(5.10) 


(5.11) 


u.tx,e>=-p«,o;-A  cf^Un'x)-^  f  j&F 


i=m‘  ?  (5  ttiY'"x>  +i  *  f  (k^ftt-3 


(5.  12) 


(5. 13) 

and  these  are  again  shown  on  Figs.  12,  13,  14  and  15.  The  conclusions  to 
be  drawn  from  these  graphs  are  broadly  similar  to  those  for  the  case  of 
the  single  wire.  Once  again  the  changes  in  velocity  are  largely  independent 
of  Mach  numbers,  but  high  Mach  numbers  again  result  in  large  expansion 
for  y  =1.4.  Note  that  in  Fig.  15a  the  density  change  for  M  =  5  is  given 
while  in  Fig.  15b  M  =  2  is  the  highest  value  plotted.  As  far  as  the  pressure 
is  concerned,  the  dominant  effect  is  again  transonic;  a  large  rise  in  pressure 
appears  asM-*l. 

Channel  flow  problems  in  the  supersonic  case  are  more  difficult 
in  view  of  the  multiple  reflections,  but  a  solution  in  the  form  of  an  infinite 
series  might  possibly  be  constructed.  No  such  cases  have  been  attempted. 

6.  Variable  Conductivity  Problems 

Apart  from  problems  involving  physical  boundaries,  an  interesting 
and  in  many  ways  realistic  problem  is  that  in  which  the  conductivity  is 
allowed  to  vary  with  position.  Up  to  this  point,  it  has  been  assumed  that  the 
conductivity  a*  was  a  constant  at  all  points.  It  is  possible  by  using  the  methods 
that  have  been  developed  in  this  paper  to  treat  problems  in  which  a*  is 
assumed  to  be  a  function  of  y  only.  Attention  will  be  restricted  to  the 
incompressible  case  with  a  single  wall  at  y  =  0;  the  extension  to  channel 
flows,  or  to  subsonic  or  supersonic  flows  is  immediate. 

It  is  necessary  first  to  modify  the  equations  used  to  describe  the 
flow  to  include  the  possibility  of  variable  conductivity.  Taking  cr*  as  a 
reference  value  of  the  conductivity,  let  a'  be  the  actual  value  as  a  function 
of  position,  and  define  a  =  a'/ a*.  Then,  Eqs.  (2.  11)  to  (2.  15)  will  be 
modified  by  replacing  S  with  oS.  Hence,  using  Eqs.  (2.  19)  to  (2.  23),  the 
right-hand  sides  of  Eqs.  (2.  24),  (2.  25),  and  (2.  27)  should  be  multiplied  by 
a.  In  particular,  the  equations  governing  p  and  v  become 
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(6.1) 


'  0  -  1^-  =  <r  D -Klf- 1 ) 


(6.2) 


For  the  incompressible  case  it  is  only  necessary  to  set  M  =  0. 


The  simplest  problem  in  which  o  is  not  constant  is  where  a  =  1  for 
0  <  y  <  H  and  a  =  0  for  y  >  H.  The  appropriate  boundary  conditions  at  y  =  H 
are  that  v  and  p  should  be  continuous.  There  will,  however,  be  a  jump  in 
u.  Let  subscript  1  refer  to  the  region  0  <  y  <  H,  and  subscript  2  refer  to  the 
region  y  >  H.  Then  the  solution  can  be  constructed  as  follows;  first,  in 
region  1,  take  the  particular  solution  given  by  Eqs,  (3.4)  and  (3.  5).  This 
may  be  referred  to  as  (vpp  Pp^).  Next,  define  the  complex  function 


(6.3) 


It  is  now  required  to  find  a  pair  of  complex  functions,  wj  (z)  and  w£  (z) 
such  that  Wj  (z)  is  analytic  in  the  region  0  <  y  <  H  while  W2  (z)  is  analytic 
in  the  entire  region  y  >  H,  and  such  that 

w^Cx*  i  H )  =  wt  <x+  i  H )  -v  Fix') 

(6.4) 

However,  it  is  convenient  to  require  that  w^  should  also  be  analytic  in  the 
region  y  <  0.  In  this  case  the  general  solution  to  this  problem  is 


(6.6) 

where  6 ^  anc*  ^  are  *arge  semicircles  having  the  line  y  =  H  for  diameter 
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and  enclosing  respectively  the  regions  y<H  and  y>H.  In  simple  cases  it 
will  be  possible  to  find  the  functions  wj  and  wg  by  inspection,  but  formulae 
(6.  5)  and  (6.  6)  are  quite  general.  At  this  stage  the  solution  in  region  1  may 
be  written  symbolically  as  (Ppj>  vpi)  +  wi  >  while  in  region  2  the  solution 
is  W£.  The  conditions  of  continuity  across  the  line  y  =  H  have  been  met, 
but  the  boundary  at  y  =  0  has  not  yet  been  introduced.  In  order  to  take 
account  of  this,  a  function  w  (z)  is  introduced,  w  (z)  is  analytic  in  the  region 
y  >0  and  satisfies: 

*0  w  ■+  +  =  6  ony*o 

(6.7) 


The  final  solution  is  now  given  by: 

^Pp,  »Vp^  0<V<H 

Wx+  w  y>H 

Single  Wire  Case 

To  illustrate  the  use  of  this  technique,  consider  first  the  case  of 
a  single  wire,  the  components  of  the  magnetic  field  being  given  by  Eq.  (3.  10). 
The  solution  (ppi  ,  vpj)  is  given  by  Eqs.  (3.  11)  and  (3.  12),  and  the  function 
F  (x)  by 

Ftx),  ^[x-(,vit<)tav\VprTr1-j  »taw'  vnrr 

»)*■ 


(6.8) 


The  functions  wj  and  w£  may  be  found  by  inspection  to  be: 

I  i  zH+l+'tZ  l 

Aim  - 


w.(z)=  i- 


_  ± 


2+1 


Z-'»(Z.H  +  t) 


wzC^)= 


-  kh 1 


2+1 


I  -  I  ^  . 


I 

4- 


Z  +  \ 


(6.9) 


The  condition  on  w  at  y  =  0  is  now  3 

I  _  I  '  I  .  .  /All!  lY1* 


(6. 10) 


UAIlAV/AJl  Vi*  YY  1*1.  y  —  v  11V  VV  1 

_  L.V.  +  I>  i^H+*  +1  _iML_ 

'  X*-M *  x^+UH  +  l 


(6.11) 
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so  that 


in  48*^ 


2HT 


(6. 12) 

completing  the  solution  of  the  problem.  The  pressure  on  the  wall  is  given 
by 


p^^kcrr  - 


x*4(*-H+i)- 


ixi„ 

x*-H 


(6. 13) 


and  is  shown  in  Fig.  17.  Clearly,  when  H  vanishes  and  the  flow  is  non¬ 
conducting  everywhere,  all  the  flow  perturbations  vanish.  However,  as 
may  be  seen  from  Fig.  17,  there  is  only  a  small  difference  between  the 
case  where  the  flow  is  conducting  only  out  to  y  =  1  and  the  case  where  it 
is  conducting  everywhere. 

The  jump  in  u  at  y  =  H  is  given  by 


which  is,  as  might  be  expected,  always  positive. 
Linear  Dipole 


(6.14) 


For  the  linear  dipole  case,  with  tJ  omponents  of  the  magnetic 
field  given  by  Eq.  (3.  33), 

c  i* »)*]+£ « '){**-  (H4  0*] 

H>0=  - [x*-KhVOT3 - 


Then: 


(6.  15) 


\  f  5z^-6iz(a.H4iH/6Hi4i0H+7) 

r  480-Uiyt[  [z-iUH+i)]3  “ 


(6. 16) 
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_  I  r  3i*+4iz(H  +  *)-(8H*+22H+l7)l 

i_  wiiulK  - (THJ3 - J 


so  that  the  condition  on  w  is : 


(6.17) 


-2«z-0  (ZH+O  f3x44.fcUHz~t)x^zH-H^fi4HSzoH»7) 

=  <5F+ip  “  WH^l  (x44  (*«♦»)*]* 


(6.  18) 


giving 


-  -  (zS3iZ-4)  .  3ZZj-fcUH  +  0‘lZ-OfcH+32>Btl7) 

W  8^?  +  *)3  +  48CH+0a'Lz+'(zh+i)J^ 


The  pressure  on  the  wall  is  thus: 


(6.  19) 


rfr  ax.-^(a^±^±i£>  +  kx++z.(|oHH8B+0x2t3(*H+0*(8H  +  ^H4Sj 
Wl  i4(xi+7)3  +  ^4CH+l)i[xa+(-l.H+l)*-]3 


(6.20) 


and  is  shown  in  Fig.  18.  Here  again  it  may  be  seen  that  the  bulk  of  the 
perturbation  to  the  pressure  on  the  wall  is  induced  very  close  to  the  wall. 

The  magnetic  field  for  this  case  falls  off  with  distance  faster  than  in  the 
single  wire  case,  and  as  a  result,  if  the  flow  is  conducting  only  out  to  about 
y  =  0.  2,  the  effect  on  the  wall  pressure  is  hardly  different  from  the  case 
in  which  the  flow  is  conducting  everywhere. 

To  conclude  this  section  it  will  be  shown  how  to  adapt  solutions  of 
the  type  given  above  to  the  case  where  the  conductivity  a  is  an  arbitrary 
function  of  y.  Define  first  p^  (x,  y,  H),  v-u  (x,  y,  H)  to  be  the  solution  to 
a  problem  with  a  given  magnetic  field  in  which  cr  =  1  for  y  <  H  and  a  =  0  for 
y  >  H.  These  are  understood  to  have  different  forms  according  to  the  sign 
of  {y  -  H).  When  H  is  infinite,  p^  and  ypj  stand  for  the  solution  to  the  prob¬ 
lem  where  the  conductivity  is  everywhere  unity.  The  solution  to  the  problem 
in  which  a  is  a  given  function  of  y  is  then: 

{p(*,V)>v(X»Y)}=  <rH  |  PHC\Y,co)?  vH(x,Y,oo)j 


(6.21) 
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To  prove  this,  note  first  that  all  the  solutions  pjj  (x,  y,  H), 
vp£  (x,  y,  H)  are  continuous  throughout  the  plane  and  give  vanishing  normal 
velocity  at  the  wall.  Hence,  the  integrated  solution,  Eq.  (6.21),  has  the 
same  properties.  It  remains  to  show  that  the  integrated  solution  satisfies 
the  appropriate  differential  equation.  Taking  Eq.  (6. 1)  for  example, 

4x-  ♦  If- *•*  <*>  {-If1  + 4^  } 

-j"  { 


(6.22) 

But 

avH  (h>v) 

d*  *  dY  0  (H<V) 

(6.23) 

Thus 

4r+  C  b*by  dMa>aH  *<r(v,t*bY 

Jy 

(6.24) 


as  required;  similarly  for  Eq.  (6.2). 
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Fig.  1  This  is  the  pressure  on  the  lower  wall  of  a  channel  of  height  H 
when  the  flow  is  impeded  by  the  magnetic  field  due  to  a  current 
flowing  in  a  single  wire  at  unit  distance  below  the  wall. 


This  is  the  pressure  gradient  on  ti 
height  H  when  the  flow  is  impeded 
current  flowing  in  a  single  wire  at 


currents  situated  unit  distam 
another.  The  strength  of  the 


6 


This  is  the  pressure  gradient  on  the  wall  when  the  flow  is 
impeded  by  the  magnetic  field  due  to  two  wires  carrying  equal 
and  opposite  currents  situated  unit  distance  below  the  wall  and 
2X  from  one  another.  The  strength  of  the  current  in  either  wire 
is  varied  with  X  in  such  a  way  that  the  field  strength  at  the  origin 
remains  fixed. 


<0 
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Fig.  5b  This  is  the  pressure  gradient  on  the  wall  when  the  flow  is  impeded 
by  the  magnetic  field  due  to  two  wires  carrying  equal  and  opposite 
currents  situated  unit  distance  below  the  wall  and  2X  from  one 
another.  The  strength  of  the  current  in  either  wire  is  varied  with 
X  in  such  a  way  that  the  field  strength  at  the  origin  remains  fixed. 


01  =/ 


This  is  the  pressure  on  a  wall  when  the  flow  is  impeded  by  the 
magnetic  field  due  to  a  current  flowing  in  a  single  wire  at  unit 
distance  below  the  wall,  v  =  1.0. 


This  is  the  pressure  on  a 
magnetic  field  due  to  a  cu 
distance  below  the  wall. 


is  the  tangential  velocity  at  the  wall  when  the  flow  is  impeded 
le  magnetic  field  due  to  a  current  flowing  in  a  single  wire  at 
distance  below  the  wall,  v  =  1.4. 


magnetic  field  due  to  a  current  flowing  in  a  single  wire  at  unit 
distance  below  the  wall. 


12a  This  is  the  pressure  on 
magnetic  field  due  to  a  ' 
wall,  y  =  1.0. 


distance  below  the  wall 


S/XM 


This  is  the  density  on  a  wall  when  the  flow  is  impeded  by  the 
magnetic  field  due  to  a  linear  dipole  at  unit  distance  below  the 


7215 


Fig.  15b  This  is  the  density  on  a  wall  when  the  flow  is  impeded  by  the 
magnetic  field  due  to  a  linear  dipole  at  unit  distance  below  the 


2 


Fig.  18 


This  is  the  pressure  on  a  wall  when  the  flow  is  impeded  by  the 
magnetic  field  due  to  a  linear  dipole  at  unit  distance  below  the 
wall  and  the  fluid  is  conducting  between  the  wall  and  the  line 


APPENDIX 


Note  Added  in  Proof 


Since  this  paper  was  compiled,  two  additional  papers  on  this  general 
subject  have  come  to  the  author's  attention;  these  are  the  works  of  Ehlers^ 
and  Morioka.  7 

Ehlers  treats  the  same  physical  problem  as  is  considered  here,  but 
considers  axially  symmetric  geometries  and  obtains  solutions  by  integral 
transform  methods.  The  results  appear  in  a  form  which  does. not  lend  itself 
to  simple  interpretation.  It  is  interesting  to  note  that  particular  solutions 
to  the  axially  symmetric  incompressible  problem  analogous  to  those  given 
here  for  two-dimensional  flow  may  be  obtained.  In  fact,  the  equations 
governing  the  axi-symmetric  case  are: 


where  r  is  the  radial  coordinate. 

3bx  .  ,  3  bv  . 

~3T-  ~  "a  * 0  >  -a”5  + 

ar  a*. 


The  magnetic  field  components  satisfy: 

T*  3r  (r 


A  particular  solution  is  then 


vp  =  2.  br 


bxdx  + 


bKbr  dx 


(Al) 


(A2) 

The  particular  example  treated  by  Ehlers  involves  the  magnetic 
field  due  to  a  circular  coil.  The  field  components  in  this  case  do  not  per¬ 
mit  the  evaluation  of  the  integrals  (Al)  and  (A2)  in  closed  form,  but  it 
seems  clear  that  simple  axisymmetric  fields  could  be  treated  in  detail 
using  the  particular  integrals  given  above. 
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Morioka,  on  the  other  hand,  treats  precisely  the  two-dimensional 
problem  treated  here,  with  the  magnetic  field  being  due  to  a  linear  dipole. 
His  method  is  somewhat  similar  to  that  of  the  present  paper,  but  the  results 
disagree;  certain  misprints  have  been  found  in  Morioka' s  paper,  but  it  is 
not  known  whether  all  the  differences  in  results  can  be  attributed  to  mis¬ 
prints  . 


-58- 


REFERENCES 


1.  Fishman,  F,  ,  Lothrop,  J.  ,  Patrick,  R.  and  Petschek,  H.  ,  "Super¬ 
sonic  Two-Dimensional  Magnetohydrodynamic  Flow,  "  Avco-Everett 
Research  Laboratory  Research  Report  39,  February  1959. 

2.  Sherman,  A.  ,  "The  Effect  of  Non-Uniform  Magnetic  Fields  on  Internal 
Flows  of  Conducting  Fluids,  "  Presented  at  the  Sixth  National  Annual 
Meeting  of  the  American  Astronautical  Society,  18-21  January  I960, 

New  York  City,  Preprint  No.  60-56. 

3.  Sherman,  A.  ,  "Viscous  Magnetohydrodynamic  Boundary  Layer,  "  The 
Physics  of  Fluids,  Vol.  4,  No.  5,  May,  1961,  pp.  552-557. 

4.  Kemp,  N.  and  Petschek,  H.  ,  "Two-Dimensional  Incompressible 
Magnetohydrodynamic  Flow  Across  an  Elliptical  Solenoid,  11  Journal 
of  Fluid  Mechanics,  Vol.  4,  Part  6,  November  1958,  p.  55 J! 

5.  Jahnke,  Dr.  Eugene  and  Emde,  Fritz,  Tables  of  Functions  with  Formulae 

.and  Curves,  Dover  Publications,  1945.  ~ 

6.  Ehlers,  F,  E.  ,  "Linearized  Magnetogasdynamic  Channel  Flow  with 
Axial  Symmetry,  "  ARS  Journal,  Vol.  31,  No.  3,  March  1961,  pp.  334- 
342. 

7.  Morioka,  S.  ,  "Two-Dimensional  Flow  of  a  Weakly  Conducting  Gas 
Past  a  Plane  Wall  Containing  a  Magnetic  Dipole,  11  J.  Phys.  Soc.  Japan, 
Vol.  16,  No.  12,  December  1961,  pp.  2544-2550. 


-59- 


UNCLASSIFIED  ,  I  UNCLASSIFIED 


•  CT\ 

£  4> 

.  a  j© 

-O^hO 

8  g?  u. 


Is* 

2  H  -a  «  jS  i  a 

v  «  w  .  „  »  55  2 

a  S  a-S  r  s  *  a 

«  .2  «  .t{  4>  >  Uh  O 

Sfe2HJ<-<U 

M*  (N  fo  hH  H  5  ^  > 


■4  \o  /-\  TJ 
2  »H  iG 

*  n  *-•  *w 

3JSS  g 
J  |B  ~a 
I  Si  a 


3  %  S  g 

5  w  ja  c  £ 
£•-«{-,  w  v 
i,  n.  k  ‘ "*  n 

/>  r'g  .  w  o 

t  o  ^  £  4J 

§  c 

3  *§b  >  £  ts  ° 

?  12  H3  q  3  4> 

* « a 
*  «j  m  * 

u  C  O  •£*  O  0 
3  O  -3  “  3  52 
3.2  e  °  c  v 

L»  Jz  S)«,H^ 
:  ^  oj  c  c  " 

3  §  E MS  g 

h‘£  ^S 

^  g  h-5ti-S 

‘  «r4.s  3-S 

2  g  .  «  ~  § 

5J  SsJ.a 

!7aii  a, "3 

3  o  _  c  a.  m 

n  $  6  8  rt  .2 


■w  o  (0  o 

§  ~5  Jg  E  J3  jj 


-s  w  S' 

f  I  ^  o 

;  §  *  -s  -  £ 

:  °  b  oj** 

QJ  3 

i  f  “  » J 

i  &-.S  -S  J 

IA  S  BJ  w  CJ 

l  «->  «£)  *rj  _ 

c  j  B  •*  S 
o  „  B  tg  ■ 

!  'C  0  '■ 

u  *j  Cu  <j  v 
!  4J  ‘3  4;  ’*3  >2 

irsii 

tr*  rt  O  ,b 

*rt  ‘3  2  w 

oj  w  tr  3 

x  u  9-o  « 

*5  vt  B  -c  - 

s  «  S  5 ° 

%  «  '  M 

V  'M  “  O.S 

2T3  U4*  i 
ft  *3  *2  O 

SaiSK 

■a«i  s 

«£■£  3  § 
4J  «  .3  O  w 

a°  u  u 

u  .S'  g  2  M 

.«  >  .3  03  O 
C  ’S  «rt  rj  W 
ffl  O  «  J  « 
c  3  0  -2  3 
•  5  T3  £.53  *3 


•  S  rt  \D 

X 

T3  "  ^  O 

a  S^u, 

.•§1^ 

“wan 

•  «  I  r/)  ui 

o-s?5®e 
.3  «  >  i*.  o 


i,  J  io 

5  «  u  "  o 
g  wx  c  S 
«,HH,2  w 

..  CL.  u*  ^ 

W  ~2  -  W  o 
*  J?  U  ^  £  o 
o  t;  cj  gc 
g  *§>  >  »  8  ° 

2  33  T3  0  3  ° 

«  w«c  w  a 

0  ir12  jT  s 

a  §  -a  “  =  = 

3  Jr  vrt  ~  ^ 

0  .2  =  o  o  »> 

U  ~  §)  41  *"*  *g 

c  ^  n  c  u  ** 
3  §  s  X-S  g 
°  _r  2  01  >.  « 
M  §  H  •S  2  '3 
i,«in.s2’§ 
B  S  .»;§ 

3  ^  jO  — <  ;« 

2  '  St  o.’rt  ; 

8  %  E  i  S'.y 

u  N  (/) 

11  O  43 

3  «  ^  2  J=  ^ 

g  S  2 » °  . . 

=  =  s  s 
*«§£•§  • 

U.  ?C  ^  C  T)  , 
3  {rt  ^  7  O  4;  ' 

2  U  g  «  O  ! 

a  ,2  £  ^  jj  . 


j=  «  « 

.  3  u  o  "o 

g  9-5  “;§• 

°  fc’Oi-'3 

»  <  "U  ! 

■p  >  w  .rt  4> 

•A  2  ^  ^ 

.3^2^ 

g-°  S5  s 

•a  £  «  « 

u  rt  04  o  u 
43  '2  aj  '2 

•iJ  c  ^  ^ 

^  2.  *1  4J 

*  *  w  g  -a 

o  w.  c  6  c 
ca  J2  ,o  w  .5 

u  (3  3  w 

■S  w  e-§  J 


rt  '«  JJ  T3  o 

a  a  JS  s  « 

a  U  "  x  S 
j»  -s  -S  3  S 
2  2  S  “  2 

fl  o  ^  cd  a 

*7?  o  u 

«i  .§■  §  a  « 
.tj  t>  .2  «  o 
c  'a  .a  >  *■* 

^3  U  M  j  U 
fl  3  O  ^2  3 
.3  T3  O.G  T3 


5r*a.S  . 

|  S  »  S 

(0  rt  (LI  — i 
4>  tfl  *** 
ij  ,  3  qj 

*!"{ 

uj:  mu 

a°ll 
•S  §ijf 

O  M  .3  P 

gl  » -6 

el 
i  §  I  i 

“  j  f  5 

a  y  >  e 


.  v  X  .► 
Ji  x  —  3 

rt  U 

: 2  ^  o  a 

M  O  ^ 

W  ,.  T3  C 
4)  W  a,  o 

S  c  2  o 
CL.  o  c 
E  ’3  j;  2 

,5  3  D.’u 

£5  §1 


c  » 

o  .5} 

fr-3  e  . 

1 1  si 

4>  S*  W  **-• 
13  .  3  «i 

»l!0x 

41  M  60  ^ 

•i!  _  'u  o 

.2  s-§  •« 

C  2  3 

§=  -I 

v  «  w 

o,M  c 

Sjs  f-J 

t  w  §a 
§g2  § 

,"»«  o 

S  g.5  e 

O  B  “  W 
>2" 
3  2  S  b 
w  a  3  .2 
-  c3  .► 
4J  rt  •-<  13 

rr  u-  .o  3 
m  o  **  *2 
»  ,.-o  s 

2  S  2  8 

CL.  o  c  .. 

6*a  ^  2 
8X  S-3 
.3  3  S.-6 

^(3  §f 


UNCLASSIFIED  .  |  UNCLASSIFIED 


